Molecular docking and machine learning affinity prediction of compounds identified upon softwood bark extraction to the main protease of the SARS-CoV-2 virus.

Molecular docking of 234 unique compounds identified in the softwood bark (W set) is presented with a focus on their inhibition potential to the main protease of the SARS-CoV-2 virus 3CLpro (6WQF). The docking results are compared with the docking results of 866 COVID19-related compounds (S set). Furthermore, machine learning (ML) prediction of docking scores of the W set is presented using the S set trained TensorFlow, XGBoost, and SchNetPack ML approaches. Docking scores are evaluated with the Autodock 4.2.6 software. Four compounds in the W set achieve a docking score below -13 kcal/mol, with (+)-lariciresinol 9'-p-coumarate (CID 11497085) achieving the best docking score (-15 kcal/mol) within the W and S sets. In addition, 50% of W set docking scores are found below -8 kcal/mol and 25% below -10 kcal/mol. Therefore, the compounds identified in the softwood bark, show potential for antiviral activity upon extraction or further derivatization. The W set molecular docking studies are validated by means of molecular dynamics (five best compounds). The solubility (Log S, ESOL) and druglikeness of the best docking compounds in S and W sets are compared to evaluate the pharmacological potential of compounds identified in softwood bark.

Machine learning prediction of 3CLpro SARS-CoV-2 docking scores.

Molecular docking results of two training sets containing 866 and 8,696 compounds were used to train three different machine learning (ML) approaches. Neural network approaches according to Keras and TensorFlow libraries and the gradient boosted decision trees approach of XGBoost were used with DScribe's Smooth Overlap of Atomic Positions molecular descriptors. In addition, neural networks using the SchNetPack library and descriptors were used. The ML performance was tested on three different sets, including compounds for future organic synthesis. The final evaluation of the ML predicted docking scores was based on the ZINC in vivo set, from which 1,200 compounds were randomly selected with respect to their size. The results obtained showed a consistent ML prediction capability of docking scores, and even though compounds with more than 60 atoms were found slightly overestimated they remain valid for a subsequent evaluation of their drug repurposing suitability.

3CLpro and PLpro affinity, a docking study to fight COVID19 based on 900 compounds from PubChem and literature. Are there new drugs to be found?

The spread of a novel coronavirus SARS-CoV-2 and a resulting COVID19 disease in late 2019 has transformed into a worldwide pandemic and has effectively brought the world to a halt. Proteases 3CLpro and PLpro, responsible for proteolysis of new virions, represent vital inhibition targets for the COVID19 treatment. Herein, we report an in silico docking study of more than 860 COVID19-related compounds from the PubChem database. Molecular dynamic simulations were carried out to validate the conformation stability of compound-ligand complexes with best docking scores. The MM-PBSA approach was employed to calculate binding free energies. The comparison with ca. 50 previously reported potential SARS-CoV-2's proteases inhibitors show a number of new compounds with excellent binding affinities. Anti-inflammatory drugs Montelukast, Ebastine and Solumedrol, the anti-migraine drug Vazegepant or the anti-MRSA pro-drug TAK-599, among many others, all show remarkable affinities to 3CLpro and with known side effects present candidates for immediate clinical trials. This study reports thorough docking scores summary of COVID19-related compounds found in the PubChem database and illustrates the asset of computational screening methods in search for possible drug-like candidates. Several yet-untested compounds show affinities on par with reported inhibitors and warrant further attention. Furthermore, the submitted work provides readers with ADME data, ZINC and PubChem IDs, as well as docking scores of all studied compounds for further comparisons.